Numerical Analysis of Coherent Many-Body Currents in a Single Atom Transistor 
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We study the dynamics of many atoms in the recently proposed Single Atom Transistor setup [A. 
Micheli, A. J. Daley, D. Jaksch, and P. Zoller, Phys. Rev. Lett. 93, 140408 (2004)] using recently 
developed numerical methods. In this setup, a localised spin 1/2 impurity is used to switch the 
transport of atoms in a ID optical lattice: in one state the impurity is transparent to probe atoms, 
but in the other acts as a single atom mirror. We calculate time-dependent currents for bosons 
passing the impurity atom, and find interesting many body effects. These include substantially 
different transport properties for bosons in the strongly interacting (Tonks) regime when compared 
with fermions, and an unexpected decrease in the current when weakly interacting probe atoms 
are initially accelerated to a non-zero mean momentum. We also provide more insight into the 
application of our numerical methods to this system, and discuss open questions about the currents 
approached by the system on long timescales. 



PACS numbers: 03.75.Lm, 42.50.-p, 03.67.Lx 



I. INTRODUCTION 

The recently proposed Single Atom Transistor (SAT) 
setup provides new opportunities to experimentally 
examine the coupling of a spin- 1/2 system with bosonic 
and fermionic modes. Such couplings form fundamental 
building blocks in several areas of physics. For example, 
atoms passing through a cavity can allow the quantum 
non-demolition (QND) readout of single-photon states in 
quantum optics and in solid state physics, such sys- 
tems occur in Single Electron Transistors Q, in studies 
of electron counting statistics Q and in the transport of 
electrons past impurities such as quantum dots 

In the SAT setup, which was motivated by the sig- 
nificant experimental advances made recently with cold 
atoms in ID a single spin-1/2 impurity atom, q, 

is used to switch the transport of a gas of cold atoms in a 
ID optical lattice (Fig. [I}. The impurity atom, which can 
encode a qubit on two internal spin states, is transparent 
to a gas of probe atoms in one spin state (the "on" state), 
but acts as a single atom mirror in the other (the "off" 
state), prohibiting transport via a quantum interference 
mechanism (Fig.QJ. Observation of probe atoms that are 
initially situated to one side of the impurity, and which 
can constitute either a ID degenerage Bose or Fermi gas, 
can then be used as a QND measurement of the qubit 
state of the impurity atom \ip q ) = a\1) + (3\ [) (see 
Fig.[TJ. 

The long coherence times associated with atoms in op- 
tical lattices allow many-body effects to contribute coher- 
ently to the transport properties over longer timescales 
than is observed in other systems where bosonic and 
fermionic modes couple to a spin 1/2 system. This pro- 
duces novel physics in which the current of atoms passing 
the impurity, especially in a regime of weak coupling be- 
tween probe atoms and impurity, is sensitive to interac- 
tions between the probe atoms [l| . These effects could be 



directly observed in experiments, for example, via mea- 
surements of the density of probe atoms on each site of 
the impurity atom as a function of time. 

In this article we present a detailed numerical analy- 
sis of these currents, making use of recently developed 
numerical methods [H| to calculate the dynamics of the 
bosonic probe atoms by directly integrating the many- 
body Schrodinger equation in ID on an adaptively trun- 
cated Hilbert space. When these currents are compared 
to analytical calculations of transmission coefficients for 
single particles passing the impurity atom and the re- 
lated currents for a non-interacting ID Fermi gas, signif- 
icant interaction effects are observed, as first discussed in 
Ref. Here we provide new insight into the time de- 
pendence of these currents, and what conclusions can be 
drawn from our numerical results on different timescales. 
We then calculate the initial currents for atoms at zero 
temperature diffusing past the impurity (where the ini- 
tial mean momentum of the ID gas, (k) t =o = 0, with 
k is the operator corresponding to the quasi-momentum 
in the lowest Bloch band and t the time), and explore 
the effects observed for different interaction strengths of 
bosonic probe atoms. We then also investigate the cur- 
rents for fermions and bosons when the probe atoms are 
initially kicked ((fe) t=0 ^ 0). This study is complemen- 
tary to the analytical study of the SAT that is given in 
recent article by Micheli et al. [TT| . 

In section [H] we discuss the basic physics of the SAT, 
and give a summary of the dynamics found in p] for 
single particles and non-interacting fermions. Then we 
present in detail the numerical techniques that we use to 
compute the exact time evolution of the many-body ID 
system. The time-dependence of the resulting currents 
is discussed in section ITTT1 followed by a presentation of 
the values of the initial steady state currents, both in the 
diffusive ((fc)t= o = 0) and kicked ({k)t— o 7^ 0) regimes. 
The conclusions are then summarised in section llVl 
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FIG. 1: A Single Atom Transistor (SAT) in a ID optical 
lattice: A single spin- 1/2 impurity atom q separately trapped 
at a particular lattice site is transparent to probe atoms 6 
in one state ("on"), but in the other acts as a single atom 
mirror ("off"). The probe atoms can either diffuse past the 
impurity site with mean initial momentum {k)t=o — or can 
be accelerated to a finite initial momentum {k)t=o / by a 
kick of strength pk- 



II. OVERVIEW 
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FIG. 2: An optical Feshbach resonance for a single spin chan- 
nel (fi = Q.O, A = A CT ): One probe atom and the impurity 
atom, in an atomic state &Q<;J|vac) which is quantised by the 
trapping potential of the lattice site, are coupled by an opti- 
cal Feshbach setup to a bound molecular state, m^vac), of 
the Born-Oppenheimer potential, V(r) (note that here the 
Born-Oppenheimer potential is modified by the trapping po- 
tential of the lattice site p^|. The coupling has the effective 
two-photon Rabi frequency £1, and detuning A. 



A. The Single Atom Transistor 



2. Single Atoms 



1. The System 



As described in section [fl we consider probe atoms b, 
which are loaded into an optical lattice [121 llH . IT3. Hfij 
with strong confinement in two dimensions, so that the 
atoms are restricted to move along a lattice in ID. The 
probe atoms are initially situated to the left of a site 
containing an impurity atom q, which is trapped inde- 
pendently (by a species or spin-dependent |16| potential), 
fixing it to a particular site while the probe atoms are free 
to move. In order to produce the "on" and "off" states 
of the SAT, we must appropriately engineer the effective 
spin-dependent interaction between the probe atoms and 
the impurity, H int = Y, a U e s,ablb qlq^- Here, b\ and gt 
are second-quantised creation operators for the b and q 
atoms respectively, obeying the standard commutation 
(anti-commutation) relations for bosons (fcrmions) and 
the site index i is chosen so that the impurity is on site 
z = 0. These interactions can be controlled using either 
a magnetic [13 El or optical Feshbach resonance. 
For simplicity we discuss the case of an optical Feshbach 
resonance, depicted in Fig. [21 Here, lasers are used to 
drive a transition from the atomic state SJgJlvac) via an 
off-resonant excited molecular state to a bound molecu- 
lar state back in the lowest electronic manifold m^vac) 
on the impurity site, i — (see Fig. The two-photon 
Rabi frequency for this process is denoted Q a and the 
Raman detuning A CT , and throughout this article we use 
units with h = 1. 



We consider initially a single probe atom passing the 
impurity. If the coupling to the molecular state is far off 
resonance (fi^ <§C |A<j|), the effect of the Feshbach reso- 
nance is to modify the interaction between the b and q 
atoms in the familiar manner, with U e g = U q b + fl^./A a . 
This can be used to screen the background interaction be- 
tween these atoms, U q t, so that the "on" state of the SAT 
(Ucs — 0) can be produced by choosing Aj = —£l^/U q b- 

If the coupling is resonant ( A^ = 0) , then the physical 
mechanism is different, and the passage of a probe atom 
b past the impurity is blocked by quantum interference. 
The mixing of the unbound atomic state and the molecu- 
lar state on the impurity site produces two dressed states 



--= (SJgjlvac) 



± m] I vac 



with energies 



£± 




(1) 



(2) 



The two resulting paths for a particle of energy e then 
destructively interfere so that when ilj 3> J, where J 
is the normal tunneling amplitude between neighbouring 
lattice sites, and U q t, = 0, the effective tunnelling ampli- 
tude past the impurity (see Fig. |3J) is 



Jeff = - 



J 2 



J 2 



0. 



(3) 
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This is reminiscent of the interference effect which un- 
derlies Electromagnetically Induced Transparency [20| . 
and corresponds to the effective interaction U B g — > oo 
required for the "off" state of the SAT. 

In Refs. pi lll|. the Lippmann-Schwinger equation is 
solved exactly for scattering from the impurity of an 
atom b with incident momentum k > in the lowest 
Bloch-band, where the energy of the particle e(k) — 
— 2Jcos(fca), with a the lattice spacing. The resulting 
transmission probabilities T(p) are in the form of Fano 
Profiles |2l|. For f2 CT ~ J these have a minimum cor- 
responding to complete reflection for e(k) = — A CT and 
complete transmission for e(k) — — A CT — ti^/U q b- For 
il a > 4:J, the transmission coefficients are approximately 
independent of k, and so complete transparency of the 
impurity atom is obtained for A CT = —fl^/U q b and com- 
plete blocking of the incident atoms for A = 0. 
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FIG. 3: The sequence as (left) a probe atom approaches the 
impurity site and is located on site i = — 1, (centre) the probe 
atom is on the impurity site, i — 0, and (right) the probe atom 
has tunnelled past the impurity and is located on site i = 1. 
Quantum interference in this process because the two dressed 
states on the impurity site, |±) = (&J<jJ|vac) ± m^jvac})/-^ 
give rise to two separate paths with equal and opposite am- 
plitude. 



3. Many Atoms 

The treatment of this system for many atoms is simi- 
lar to the single atom case, but the motion of the probe 
atoms in the lattice, except on the impurity site, is gov- 
erned by a (Bose-) Hubbard Hamiltonian [F|]. As the two 
spin channels for the impurity atom, q can be treated in- 
dependently, we will consider only a single spin channel 
q a , and drop the subscript in the notation throughout 
the remainder of the article Q . The Hamiltonian for the 
system is then given (with h = I) by H = H b + Hq, with 

(y) J 
Hq = Q (rh^qbo + h.c^ — Am'tii 

+U qb blq s qb + U bm blm^rhb . (4) 

Here, H^ gives a Hubbard Hamiltonian for the b atoms 
with tunnelling matrix elements J, and collisional in- 
teractions Ubb- For fermions, XJbb — 0, whereas for 
bosons Ubb = 4:Trh 2 abb J d 3 x. |wj(x)| /rrib, with Wj(x) the 
Wannicr- function on site j, and abb and rrib the scatter- 
ing length and mass of b atoms respectively. i?o describes 
the dynamics in the presence of the impurity on site 0, 
where atoms b and q are converted to a molecular state 
with effective Rabi frequency and detuning A, and the 
final two terms describe background interactions, U a js for 
two particles a, (3 € {q,b,m}, which are typically weak 
and will be neglected in our treatment. This single-band 
model is valid in the limit for U a p, J, fi,A <C w, where 
lu is the energy separation between Bloch bands, an in- 
equality which is fulfilled in current experiments. The 
robustness of the SAT with respect to loss processes is 
discussed in . 

In the rest of this article, we will study the current of 
atoms past the impurity site that develops as a function 
of time, and how this current depends on the interaction 



between probe atoms and on interactions between the 
probe atoms and the impurity. 



B. Atomic Currents through the SAT 

To analyse the case of many atoms passing the impu- 
rity site, we consider the probe atoms b to be prepared 
initially to the left of the impurity, in a ground state 
corresponding to a ID box potential. The current of 
atoms passing the impurity is I(t) — dNn/dt, where Nr 
is the mean number of atoms to the right of the impurity, 
^fi = (J2j>ob]bj)- For a sufficiently large number of 
atoms in the initial cloud, this current is generally found 
to rapidly settle into an initial steady state current, I ss , 
on relatively short timescales (tJ ~ 1) (see section IlII Al 
for further discussion of steady state currents for bosons). 

For a non-interacting Fermi gas at zero temperature, 
the currents can be calculated exactly when Ubm — U q b, 
as the equations of motion are linear. Scattering from 
the impurity occurs independently for each particle in 
the initial Fermi sea, and after a short transient period 
of the order of the inverse tunnelling rate 1/J, a steady 
state current I ss is established. This can be calculated ei- 
ther by int egra ting the single-particle transmission prob- 
abilities 0, El or by direct numerical integration of the 
Heisenberg equations. 

For a non-interacting and very dilute Bose gas, the 
situation will be identical to considering a single parti- 
cle. However, for higher densities, many-boson effects 
become important, and additionally for non-zero inter- 
actions the situation becomes even more complicated. In 
the limit Ubb/ J — > oo in fD (the Tonks gas regime) it is 
usually possible to replace the bosonic operators bi , b \ by 
fermionic operators /»,/,• using a Jordan- Wigner trans- 
formation 22]. However, in this case the resulting Hamil- 
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tonian, 

m 

+Uqbfotfqfo + Ubmfom^rhfo, (5) 

contains a nonlinear phase factor resulting from the 
coupling on the impurity site, (— 1) L , where Nl — 

^j<o fj fi ^ s ^ ne °P era tor for the number of atoms to the 
left of the impurity site. For £1 = 0, the boson currents 
are exactly the same as the currents for noninteracting 
fermions as (Sj&i) = (/}/%)■ For finite f2 it is not clear 
what role the phase factor will play in determining the 
system dynamics, although for sufficiently large Q 3> J 
we again expect very little current to pass the impurity. 

Thus, for the intermediate regime f2 ~ J, and for the 
case of finite interaction strength U / J there are no known 
analytical solutions for the currents. For this reason, we 
specifically study these regimes in this paper, using near- 
exact numerical methods. 



C. Time-Dependent Numerical Algorithm for ID 
Many-Body Systems 

The algorithm that we use to compute the time evo- 
lution of our many body system for bosonic probe 
atoms was originally proposed by Vidal [lOj . This 
method allows near-exact integration of the many body 
Schrodinger equation in ID by an adaptive decimation of 
the Hilbert space, provided that the Hamiltonian couples 
nearest-neighbour sites only and that the resulting states 
are only "slightly entangled" (this will be explained in 
more detail below). Recently both this algorithm [23| . 
and similar methods proposed by Verstrate and Cirac 
[24j have been generalised to the treatment of master 
equations for dissipative systems and systems at finite 
temperature, and progress has been made applying the 
latter method to 2D systems plf . 

In ID, these methods rely on a decomposition of the 
many-body wavefunction into a matrix product repre- 
sentation of the type used in Density Matrix Renormali- 
sation Group (DMRG) calculations [2(j , which had pre- 
viously been widely applied to find the ground state in 
ID systems. The time dependent algorithms have now 
been incorporated within DMRG codes |27(, and also 
been used to study the coherent dynamics of a variety 
of systems @. In our case, we write the coefficients 
of the wavefunction expanded in terms of local Hilbert 
spaces of dimension S 1 , 

s 

I*) = Yl Ciiia...wl*l>®l*2)®-"®l*Af>. (6) 



as a product of tensors 

x 

r ... = V rW^iWrt 2 ]* 2 A Mr M* 2 H M l iM 

ai ...a M -i 

(7) 

These are chosen so that the tensor A« specifies the coef- 
ficients of the Schmidt decomposition (23 for the bipar- 
tite splitting of the system at site I, 

m = f:^-- l wr--- M] ), (8) 

a=l 

where xi is the Schmidt rank, and the sum over remain- 
ing tensors specify the Schmidt eigenstates, |</>£"'^) and 
\(f)a +1 "' M ^) . The key to the method is two- fold. Firstly, 
for many states corresponding to a low-energy in ID sys- 
tems wc find that the Schmidt coefficients A«, ordered 
in decreasing magnitude, decay rapidly as a function of 
their index a (this is what we mean by the state being 
"slightly entangled") ^(j- Thus the representation can 
be truncated at relatively small \ and still provide an 
inner product of almost unity with the exact state of the 
system \^>). Secondly, when an operator acts on the local 
Hilbert state of two neighbouring sites, the representa- 
tion can be efficiently updated by changing the T tensors 
corresponding to those two sites, a number of operations 
that scales as x 3 ^ 3 f° r sufficiently large x JO). Thus, we 
represent the state on a systematically truncated Hilbert 
space, which changes adaptively as we perform opera- 
tions on the state. 

In order to simulate the time evolution of a state, 
we perform a Suzuki- Trotter decomposition 30] of the 
time evolution operator exp(— iHi), which is applied to 
each pair of sites individually in small timesteps St. Ini- 
tial states can also be found using an imaginary time 
evolution, i.e., the repeated application of the operator 
exp(— HSt), together with renormalisation of the state. 

In this paper, results are not on ly p roduced using the 
original algorithm as presented in [HJ], but also using an 
optimised version in which the Schmidt eigenstates are 
forced to correspond to fixed numbers of particles. This 
allows us to make use of the total number conservation 
in the Hamiltonian to substantially increase the speed 
of the code, and also improve the scaling with x an d 
S. With this number conserving code we are able to 
compute results with much higher values of however 
we also find that for insufficiently large Xi the results 
from this code become rapidly unphysical, in contrast to 
the original code (see section UlI All . 

In implementations of this method we vary the value of 
X to check that the point at which the representation is 
being truncated does not affect the final results. A useful 
indicator for convergence of the method is the sum of the 
Schmidt coefficients discarded in each time step, although 
in practice the convergence of calculated quantities (such 
as the single particle density matrix, (b\bj)) are normally 
used. This is also discussed further in section ITlI AI 
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For bosons on an optical lattice we must also choose 
the dimension S of the local Hilbert space, which corre- 
sponds to one more than the maximum number of atoms 
allowed on one lattice site. For simulation of the SAT, we 
allow a variable dimension of the local Hilbert space Si , 
as we must consider the state of the molecule on the im- 
purity site in addition to the probe atoms. Allowing such 
a variable dimension dramatically reduces the simulation 
time, which scales as x 3 Yli when x ^> S, and scales 
proportional to S 4 when x is small. For a Bose gas with 
finite U/J we usually take Si — 6 away from the impurity 
site, and So = 12 on the impurity site, whereas simula- 
tions of a Tonks gas can be performed with Si =2 away 
from the impurity site and So = 4 on the impurity site. 



III. NUMERICAL RESULTS 

In section ITU Al we discuss the time dependence of the 
current for bosons and the applicability of our numerical 
methods in different regimes. We establish the existence 
of an initial steady state current, Iss that appears on a 
timescale tJ ~ 1, and discuss the observation of a sec- 
ond steady state current Iq, observed in some cases on 
a timescale tJ ~ 10. In sections IlII Bl and IIII CI we then 
present our numerical results for Iss for the case where 
the initial cloud diffuses past the impurity site, and the 
case where the initial cloud is kicked respectively. 

We are primarily interested in the behaviour of the 
current through the SAT when it is used in the "off" 
state, i.e., we choose A = 0. To enhance clarity of the 
results, we also choose Ub q = Ubm = 0. 

In each case, we considered an initial cloud of between 
N = 1 and N = 30 atoms, confined on M = 30 lattice 
sites situated immediately to the left of the impurity site. 
The initial state used corresponds to the ground state, 
\<po) of a Bose- Hubbard model with a box trap. 

Our total grid for the time evolution consisted of 61 
lattice sites, with the 30 rightmost sites initially unoccu- 
pied, and the results we present are, except for very small 
systems, independent of the size of the initial cloud and of 
the grid size. Fermionic results are derived from exact in- 
tegration of the Heisenberg equations of motion, whereas 
bosonic results are near-exact simulations as described in 
section Hi CI 



A. Time Dependence of the current for bosonic 
probe atoms 

The mean number of probe atoms on the right of the 
impurity, Nr is plotted as a function of time, t, in Fig. 01 
for a Tonks gas (U/J — > oo) with fi/J = 1 and initial 
state of density n = N/M — 1. These results were cal- 
culated with the original simulation algorithm, and it is 
clear from the figure that the current settles into an initial 
steady state value Iss on the timescale t J ~ 1. However, 
as is typical for bosonic probe atoms with n > 0.5, there 
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tJ 



FIG. 4: The number of atoms to the right of the impu- 
rity site, Nr as a function of dimensionless time tJ for a 
Bose gas in the Tonks limit (Ubb/J — > oo) with Q./J = 1, 
n = N/M — 1, and varying number of states retained in the 
method, x = 10,20,30,40,50,60,70 (lines from bottom to 
top). These results are from the original simulation method. 



exists a knee in the curve at a time £knee(x)j leading to 
a new and final steady state current, which we will de- 
note Jo- The time iknec(x) depends on the initial density, 
n, and coupling, f2, and as can be seen from this figure, 
we require a high value of x to find the exact time. For 
n = 1,0/ J = 1, iknce(x) appears to converge to a value 
between t J = 9 and t J = 12 as x is increased. It is clear 
that significant level of correlation, or entanglement be- 
tween the left and right hand side of the system (in the 
sense of the number of significant Schmidt eigenvalues for 
a bipartite splitting) are involved in determining the dy- 
namics leading the to knee. However, the actual value of 
the steady state current Iq appears to converge for much 
lower values of x and there is essentially no change in 
this result from x = 10 to \ — 70. 

The interpretation of these results is more complex 
when they are compared with similar results from the 
new, number conserving version of our code. In Fig. [5] 
we observe that the behaviour diverges at the same value 
°f ^knoc(x); and even for x = 300, the value of ikncc(x) 
has only shifted a little further from where it was ob- 
served for x — 70 with the original version of the code. 
This confirms that the dynamics on this timescale are 
dominated by the significant level of correlation, or en- 
tanglement between the left and right hand side of the 
system. 

In contrast to the steady state current Iq obtained us- 
ing the original code, though, the current in the num- 
ber conserving simulations rapidly approaches 0, even 
for x = 300. As can be seen from the dotted line in 
Fig. [BJ this behaviour occurs when the maximum sum 
of squares of the Schmidt coefficients being discarded in 
each timestep, e\ = J2p> x ^% reaches a steady value on 
the order of 10 -7 , indicating that the simulation results 
from the number conserving code are probably not valid 
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FIG. 5: The number of atoms to the right of the impurity site, 
Nr as a function of dimensionless time t J for a Bose gas in the 
Tonks limit (U bb /J -> oo) with SI/ J = l,n = N/M = 1. This 
plot shows a comparison of results from the original method 
(dashed lines, \ — 50,70, c.f. Fig. and from the number 
conserving method (solid lines, \ — 50, 100, 200, 300). 




FIG. 6: Comparison as a function of dimensionless time tJ 
of the number of atoms to the right of the impurity site, Nr, 
(solid line), and the sum of squares of the discarded Schmidt 
eigenvalues, e\ = J2/3 >x ^3 (dashed line). These results are 
taken from the number conserving simulation method with 
X = 100, for a Bose gas in the Tonks limit (Ubb/J ~* oo) with 
SI/ J = 1, n = N/M = 1. 



for t > ikncc- Indeed, we observe the same behaviour 
from the new simulation code with = 0, where we 
know from Eq. that the time dependent current I(t) 
is equal to that for fermions, and should not decrease in 
this manner (see currents for fermions in Ref. In- 
terestingly, the original code, which produces the steady 
state currents Iq at finite SI reproduces the known result 
at SI = exactly even for small values of %, with a steady 
state current Iss and no knee. 

Our conclusions from these results are as follows: 

(i) We know that up to iknoc our simulation results are 
exact, as they are unchanged in the linear region with 
current Iss for \ = 20 — > 300. As this regime lasts at 
least until t J ~ 10, these results would be observable in 
an experimental implementation of the SAT. 

(ii) As an impractically large value of x would 
be required to reproduce the results exactly on long 
timescales, we can not be certain what the final behaviour 
will be for t > ikncc(x — 300). This depends on clearly 
interesting phenomena that arise from strong correlations 
between the left and right sides of the impurity site, and 
could include settling to a final steady state current Iq. 
These effects would also be observable in an experiment. 

The expected final steady state values Jo are already 
discussed in Ref. and so in the remainder of this 
article we investigate the initial steady state currents Iss 
in various parameter regimes. 



B. Diffusive evolution, with initial mean 
momentum ((fc) t = = 0) 

We first consider the motion of atoms past the impurity 
site in the diffusive regime, where the initial state at t = 



is the ground state of a Bose-Hubbard model on M = 30 
lattice sites in a box trap. 



1. Dependence of the current on impurity-probe coupling, SI 

In Fig. [7] we show the initial steady state current Iss 
as a function of SI/ J for fermionic probe atoms, and for 
bosonic probe atoms with Ubb/J = 4, 10, oo and A = 0. 
All of these results decrease as expected with increasing 
SI/ J, and even for a relatively small SI = 2 J the current is 
minimal in each case. At half filling (Fig. [7^,), the results 
for the Tonks gas are identical to the Fermi results for 
SI = 0, but become substantially different as SI increases, 
with the currents in this regime greater for the bosons. 
At weaker interactions the currents are smaller than the 
Tonks result at all SI, but for SI/ J > 1 the currents for 
U/J = 4 are larger than for a non-interacting Fermi 
gas. The variation in the currents for different interac- 
tion strengths of bosons appears to be due to the broader 
initial momentum distributions that occur at larger U/J. 
At unit filling (Fig. UJ)), Iss is less dependent on the in- 
teraction strength, with all of the bosonic results very 
close to one another, currents becoming larger than that 
for fermions when SI/ J > 1 . 



2. Dependence of the current on interaction strength, U/J 

The dependence of the initial steady state current Iss 
on the interaction strength for bosons is depicted more 
clearly in Fig. [5] both at unit filling, n = 1, and half 
filling, n — 1/2 for SI = J. At half filling the cur- 
rent increases with increasing U/J, which is due to the 
broader initial momentum distribution produced by the 
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FIG. 7: Steady state currents through the SAT Iss as a 
function of the coupling between probe atoms and the im- 
purity, Q./.J. These plots show the comparison of a Bose gas 
with different interaction strengths U/J — 4 (dotted line), 
U/J — 10 (dashed) and U/ J — > oo (solid), and a Fermi gas 
(dash-dot), with (a) n = 1/2 and (b) n = 1. In both cases, 

A = t/qi, = (7 bm = 0. 




FIG. 8: Steady state currents through the SAT I aa as a func- 
tion of the interaction strength U /J for bosonic probe atoms 
initially at half- filling, n = 1/2 (squares) and at unit filling, 
n — 1 (diamonds) with S7/J = 1. The equivalent results for 
the Tonks gas (U/J — ► oo) and fermions are marked on the 
right hand site of the plot. A = U q b = Ubm = 0. 



higher interaction energies. In contrast, at higher den- 
sities (here n = 1), the probe atoms are blocked better 
by the SAT for higher interaction strengths, and Iss de- 
creases. The key principle here is that bosons appear to 
be better blocked when they approach the impurity in- 
dividually. For high densities this is achieved when large 
interaction strengths eliminate the higher occupancies of 
all lattice site including the impurity site. For weaker 
interactions the bosons can swamp the transistor, with 
one atom being bound to the impurity, whilst other probe 
atoms tunnel onto and past the impurity site. 

This effect is seen in Fig. EI where the molecular occu- 
pation and average probe atom occupation on the impu- 
rity site are shown for (a) U / J = 4 and (b) U/J — 10. We 
see that as n increases, the molecular occupation becomes 
rapidly higher for U/J = 10 than for U/J = 4, despite 
the larger occupation of probe atoms on the impurity site 
for U/J = 4. This indicates that for U/J = 10 atoms 
arrive individually at the impurity site, where they are 
coupled with the impurity atom into a molecular state, 
and their transport is efficiently blocked. For U/J = 4, 
more than one atom enters the impurity site at once, 



n„ 0.3 




FIG. 9: Plot showing the average steady state occupation 
of the molecular state (dashed lines) and the average steady 
state atomic occupation of the impurity site (solid lines) for 
(a) U/J = 4 and (b) U/J = 10, as a function of (a, b 1) 
0,/J with n — 1 and (a, b 2) n with f2/J = 0.5. In all cases, 
A = Uqb = Ubm = 0, and calculations were performed for 
M = 30. 



leading to a larger average probe atom occupation on 
the impurity site, but a comparatively small molecular 
occupation. 

It is important to note, however, that even when 
U/J = 4, the resulting currents are only slightly larger 
than they are for non-interacting fermions. At higher 
interaction strengths we then see an even stronger sup- 
pression of the steady state current for dense, strongly 
interacting bosons. As f2 increases, both the molecular 
occupation and probe atom occupation on the impurity 
site decrease (Fig. |5J) as the probability of even a single 
atom tunnelling onto the impurity site becomes small. 
For £1 > 2J the blocking mechanism of the SAT func- 
tions extremely well even in the regime where the probe 
atoms are dense and weakly interacting. 



3. Dependence of the current on initial density, n 

In Fig.^Jwe show the dependence of the initial steady 
state current, Iss on the initial filling factor n with (a) 
Q.J J = 0.5 and (b) Q/J = 1. In both cases, the currents 
for bosons of different interaction strengths are very simi- 
lar, with the variations following the patterns discussed in 
the preceding section. These results also agree well with 
the results for fermions at small n and for n ~ 1, but the 
plateau observed in fermionic currents near n ~ 0.5 docs 
not occur in the currents for bosons. For fermions, this 
plateau arises from the transmission profile of the SAT as 
a function of incoming momentum [lj , and occurs when 
the Fermi momentum is raised past the minimum in this 
transmission profile. For interacting bosons, this corre- 
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(b) 















0.2 0.6 1 0.2 0.6 1 
n n 

FIG. 10: Steady state currents through the SAT I sa as a 
function of the initial density of atoms n — N/M. These plots 
show the comparison of a Bose gas with different interaction 
strengths U/J = 4 (dotted line), U/J = 10 (dashed) and 
U/J — > oo (solid), and a Fermi gas (dash-dot), with (a) Q/J = 
0.5 and (b) Q/J = 1. In both cases, A = U q t = Ubm = 0. 




FIG. 11: Steady state currents I$s through the SAT for 
fermions as a function of the kick parameter pt, for (a) Q — 
and (b) Q/J=l. In each plot the lines from bottom to top se- 
quentially correspond to filling of N = 3, 6, 9, 12, 15, 20, 25, 30 
particles initially on M — 30 lattice sites. Note that the scales 
are different for (a) and (b), and also that the results in (a) 
are exactly the same as those for a Tonks gas of bosons. In 
both cases, A = U q b = Ubm = 0. 



spondence between the momentum distribution of the 
gas and the transmission profile is destroyed by many- 
body effects, and we see instead a smooth increase in the 
current. This results in the bosonic currents being sub- 
stantially larger than those for fermions near half filling 
when f2 ^ 1 (as was previously observed in Fig. Up,). 

C. Kicked evolution, with initial mean momentum 

{{k}t=o ± 0) 

In this section we consider an initial state with a non- 
zero initial momentum, which is obtained, e.g., by briefly 
tilting the lattice on a timescale much shorter than that 
corresponding to dynamics of atoms in the lattice. If the 
tilt is linear, the resulting state will be given by 

W = 0)> = ^expfafcjStSjOlfo), (9) 

3 

where |</>o) is the initial many-body ground state, and 
the quantity pk is determined by the magnitude and du- 
ration of the tilt. The effect of this tilt is to translate 
the ground state in the periodic quasimomentum space 
by a momentum pk- The final mean momentum (fc) then 
depends both on the value pk and the properties of the 
initial momentum distribution. 



1. Dependence of the current on kick strength pk 

In the case of fermions, the dependence of the cur- 
rent on q for different filling factors n — N/M and f2 
can be clearly understood in terms of the SAT transmis- 
sion profile (see ^2). ^ n Fig. II lb we see the current I ss 
as a function of pk with O = 0. The currents are each 
peaked at pk — 7r/2, where the resulting mean velocity 
of the probe atoms is the largest. For N/M = 1, the 
whole Bloch band is filled, and the momentum distribu- 
tion is not changed by the application of the kick, i.e., 



(k)t=o = 0. In Fig. Illb the same results are shown, but 
with Q/J = 1. Here we see that for small filling factors, 
a minimum appears at pk = tt/2, corresponding to the 
minimum in the transmission profile of the SAT for this 
incident momentum 0,0]. At higher filling factors, this 
feature of the transmission profile for 0/ J = 1 is not suf- 
ficiently broad to overcome the increase current due to 
higer mean velocities in the initial cloud, and the peak 
at pk — 7r/2 reappears. The currents here are, of course, 
reduced in comparison with those for O = 0. 

Whilst for all pk the currents with no coupling to the 
impurity atom, i.e., O = 0, are the same for the Tonks 
gas as for fermions (Fig. Ill^ t). the currents for finite in- 
teraction strengths are found to be remarkably different. 
In Fig. El these rates are plotted for U/J — 4,7, 10 for 
N = 5, 15, 30 particles initially situated on M = 30 sites. 
For the very dilute system with N = 5 (Fig. 112k ) we 
see a peak similar to that observed for fermions which 
is independent of the interaction strength. Here the cur- 
rents are essentially those for non-interacting particles, 
and the currents determined by the initial momentum 
distribution. For N = 15 (Fig. 112b) we observe the sur- 
prising result that the current is peaked at a lower value 
than is observed for fermions or for the Tonks gas. We 
have observed this peak consistently for such cases of fi- 
nite interaction, and note that as U / J increases, the peak 
moves back towards pk — 7r/2 as the currents converge to 
the Tonks gas results. As N is further increased, the peak 
continues to move left, and for TV = 30 (Fig. [TSb l we see 
a monotonically decreasing current as pk increases. As 
U/J increases these values tend towards the pk indepen- 
dent result observed for the Tonks gas. These results are 
surprising, but the trends in the behaviour are clear, and 
they should be directly verifiable in experiments, even 
without the presence of the impurity atom. 

For non-zero coupling to the impurity atom, the cur- 
rents as a function of pk are shown in Fig. 1131 Again we 
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FIG. 12: Steady state currents with coupling to the SAT, 
Q, = 0, Iss as a function of the kick parameter pt for varying 
interaction strengths, U/J = 4 (dotted), U/J — 7 (solid), 
and U/J = 10 (dashed), for (a) TV = 5, (b) TV = 15 and (c) 
TV = 30 particles initially situated on M — 30 lattice sites. In 
all cases, A = U q b = Ubm = 0. 




FIG. 13: Steady state currents through the SAT Iss as a 
function of the kick parameter p^, for fl/J — 1. These plots 
show the comparison of a Bose gas with different interaction 
strengths U/J = 4 (dotted line), U/J = 10 (dashed) and 
U/J — ► oo (solid), and a Fermi gas (dash-dot), with (a) n = 
1/2 and (b) n = 1. In both cases, A = U q b = Ubm = 0. 



notice that the current for bosons with finite interaction 
strength is peaked at much lower values of pk than the 
fcrmionic currents and that peaks of all of the bosonic 
currents, including the Tonks currents, as significantly 
larger than the fermionic currents at half filling, as was 
observed for diffusive results (pk = 0). The most remark- 
able feature of these plots is that despite a significant 
reduction in the current, the basic dependence on p^ is 
very similar to the O = results. 



0.7 



0.35 




FIG. 14: A comparison of steady state currents through the 
SAT, Iss, as a function of fl/J for p^ = (solid line), n/A 
(dashed), and n/2 (dotted). Here we consider a Tonks gas 
(U/J — * oo) of bosonic probe atoms which is initially at half 
half-filling, TV = 15, M = 30. A = U qb = U bm = 0. 



2. Dependence of the current on impurity-probe coupling, 

The steady state current I ss is shown in Fig. [21 as a 
function of f2. We observe the same strong decrease in 
the current due to the operation of the SAT for all of 
these curves, with the highest currents corresponding to 
the pk = tt/2 curve as expected. Note that the value of 
Iss is affected equally for all of the kick strengths, and 
the ratio in the currents for different values of pk is very 
similar for £1 = and 0/J= 2. 



IV. SUMMARY 

In summary, the SAT setup provides new experimen- 
tal opportunities to study coherent transport of many 
atoms past a spin-1/2 impurity due to the relatively long 
coherence times that exist for systems of atoms in op- 
tical lattices. The resulting coherent many-body effects 
can be clearly seen in the difference between the atomic 
currents observed for fermions and bosons, and the non- 
trivial dependence of the current on interaction strength 
for bosons with finite interactions. Even stronger depen- 
dence on these interactions is observed when the probe 
atoms are initially accelerated to a non-zero momentum. 
The initial steady state currents would be directly ac- 
cessible quantities in the experimental implementation of 
the SAT, and using recently developed methods for time- 
dependent calculation of many-body ID systems, we have 
made quantitative predictions for the corresponding cur- 
rents for a wide range of system parameters. We cannot 
be certain about the values the currents approach at long 
times, although it is possible that the system will set- 
tle eventually into a regime with a different steady state 
current. The high values of x needed to reproduce this 
behaviour in our numerical calculations suggest that the 
currents in this regime could also be strongly sensitive to 
the coherence properties of the system, which would be 
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very interesting to investigate in an experiment. 
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